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ABSTRACT 

Context. The internal dynamics of ultra-compact dwarf galaxies (UCDs) has attracted increasing attention, with most of the UCDs 
studied to date located in the Virgo cluster. 

Aims. Our aim is to perform a comprehensive census of the internal dynamics of UCDs in the Fornax cluster, and to shed light on the 
nature of the interface between star clusters and galaxies. 

Methods. We obtained high-resolution spectra of 23 Fornax UCDs with -10.4 > M v > -13.5 mag (10 6 < M/Mq < 10 s ), using 
FLAMES/Giraffe at the VLT. This is the largest homogeneous data set of UCD internal dynamics assembled to date. We derive 
dynamical M/L ratios for 15 UCDs covered by HST imaging. 

Results. In the M v -cr plane, UCDs with M v < -12 mag are consistent with the extrapolated Faber- Jackson relation for luminous 
elliptical galaxies, while most of the fainter UCDs are closer to the extrapolated globular cluster (GC) relation. At a given metallicity, 
Fornax UCDs have, on average, M/L ratios lower by 30-40% than Virgo UCDs, suggesting possible differences in age or dark matter 
content between Fornax and Virgo UCDs. For our sample of Fornax UCDs we find no significant correlation between M/L ratio 
and mass. We combine our data with available M/L ratio measurements of compact stellar systems with 10 4 < M/Mq < 10 8 M, 
and normalise all M/L estimates to solar metallicity. We find that UCDs (M>2x10 6 Mq) have M/L ratios twice as large as GCs 
(M<2x10 s Mq). We argue that dynamical evolution has probably had only a small effect on the current M/L ratios of objects in the 
combined sample, implying that stellar population models tend to under-predict dynamical M/L ratios of UCDs and over-predict those 
of GCs. Considering the scaling relations of stellar spheroids, we find that UCDs align well along the 'Fundamental Manifold'. UCDs 
can be considered the small-scale end of the galaxy sequence in this context. The alignment for UCDs is especially clear for r e > 7pc, 
which corresponds to dynamical relaxation times that exceed a Hubble time. In contrast, globular clusters exhibit a broader scatter 
and do not appear to align along the manifold. 

Conclusions. We argue that UCDs are the smallest dynamically un-relaxed stellar systems, with M>2x10 6 Mq and 7<r c /pc<100. 
Future studies should aim at explaining the elevated M/L ratios of UCDs and the environmental dependence of their properties. 

Key words, galaxies: clusters: individual: Fornax - galaxies: dwarf - galaxies: fundamental parameters - galaxies: nuclei - galaxies: 
star clusters 



1. Introduction 

In recent years, significant effort has been devoted to studying 
the internal dynamics of extragalactic compact stellar systems in 
the mass regime of massive globular clusters and ultra-compact 
dwarf galaxies (10 6 < M/M Q < 10 8 ) (Drinkwater et al. 2003, 
Martini & Ho 2004, Hasegan et al. 2005, Maraston et al. 2004, 
Rejkuba et al. 2007, Evstigneeva et al. 2007, Hilker et al. 2007). 
A compilation of the available data is presented in Mieske & 
Kroupa (2008) and Dabringhausen, Hilker & Kroupa (2008). 

A striking outcome of these studies is the finding that the 
dynamical M/L ratios of massive compact stellar systems are, 
on average, about two times larger than those of normal globu- 
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lar clusters of comparable metallicity. Several objects have M/L 
ratios at the limit of, or even beyond, the range predicted by 
stellar population models assuming canonical IMFs (Hasegan et 
al. 2005). Possible explanations for these high M/L ratios in- 
clude extreme stellar mass functions (Mieske & Kroupa 2008; 
Dabringhausen, Hilker & Kroupa 2008) or densely packed dark 
matter (Goerdt et al. 2008). The occurence of objects with high 
M/L ratios is observed to start at ~ 2x10 6 Mq (Hasegan et al. 
2005, Rejkuba et al. 2007, Mieske & Kroupa 2008), coinciding 
with "breaks" in physical size (Hasegan et al. 2005, Mieske et 
al. 2006, Kissler-Patig et al. 2006) and stellar content (Mieske et 
al. 2006). These breaks are consistent with the hypothesis that 
at ~ 2xlO 6 M0 (My — — 11 mag) we observe the transition be- 
tween simple globular clusters and more complex systems, the 
UCDs. 
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Three main formation scenarios have been suggested for 
UCDs: 1. UCDs are stellar super clusters formed in the tidal 
arms of violent gas-rich galaxy mergers (Fellhauer & Kroupa 
2002, 2005). 2. UCDs are tidally stripped compact remnants of 
nucleated dwarf galaxies (Bassino et al. 1994, Hiiker et al. 1999, 
Bekki et al. 2003, Goerdt et al. 2008). 3. UCDs are genuine com- 
pact dwarf galaxies formed from small-scale peaks in the pri- 
mordial dark matter power spectrum (Drinkwater et al. 2004). In 
the first case, UCDs are not expected to contain any dark matter. 
In the second and third case, UCDs would be related to cosmo- 
logical low-mass dark matter halos and may contain dark matter 
(Goerdt et al. 2008). The comparably large M/L ratios of 5-10 
found for some Virgo UCDs may point towards a cosmological 
origin. These M/L ratios are similar to the values found for some 
of the more luminous Local Group dSphs like Sculptor and Leol 
(Gilmore et al. 2007), although note that they are still 1-2 or- 
ders of magnitude below the M/L values found for the ultra-faint 
dSph candidates (Gilmore et al. 2007, Simon & Geha 2007). 

In this paper, we aim to study whether the high M/L ratios are 
a fundamental trend equally common to all UCDs, or whether 
environmental variations of the trend exist. The latter may be ex- 
pected if competing formation channels dominate in different en- 
vironments (e.g. Mieske et al. 2006). For both the Virgo cluster 
(Hasegan et al. 2005, Evstigneeva et al. 2007) and the Centaurus 
A group (Rejkuba et al. 2007), more than 10 objects in the UCD 
mass range have measured M/L ratios. For the Fornax cluster, 
only five sources have high resolution spectroscopy available 
(Hilker et al. 2007), making it difficult to judge whether differ- 
ences exist between Fornax and Virgo UCDs. With a compre- 
hensive sample of measurements for Fornax, it will be possible 
to analyse whether Fornax UCDs extend to such high M/L ratios 
as Virgo UCDs. 

In what follows, we present new measurements of the inter- 
nal kinematics of 23 compact objects in Fornax (obtained with 
FLAMES/GIRAFFE at the VLT). We analyse how their M/L 
ratios relate to predictions from stellar population models, and 
investigate how they fit into the trend of increasing M/L with 
mass among compact stellar systems. We also examine how 
UCDs and globular clusters fit into the broader context of larger 
and more luminous stellar systems, focusing on the fundamen- 
tal manifold of stellar spheroids (Zaritsky et al. 2006a, 2006b, 
2008) that is an extension of the fundamental plane concept. 

The paper is structured as follows: in Sect. [2] we present the 
new spectroscopic data of Fornax UCDs used for the present 
study. Section|3]describes the data reduction, including the mod- 
elling of the mass distribution. In Sect. @] the results for the 
Fornax UCDs are presented and discussed. In Sect. [5] we com- 
bine the Fornax data with other literature results on dynamical 
M/L ratios of compact stellar systems, and investigate UCDs and 
GCs in the context of the fundamental manifold. The paper fin- 
ishes with Summary and Conclusions in Sect. [6] Throughout this 
paper we assume a distance modulus to Fornax of (m-M)=31.4 
mag (Freedman et al. 2001). 

2. Data 

The data for this study were obtained in service mode with 
the Fibre Large Array Multi Element Spectrograph (FLAMES; 
Pasquini et al. 2002) mounted on UT2 at the VLT (programme 
078.B-0496). We used the spectrograph GIRAFFE in MEDUSA 
mode, which allows the observation of up to 130 targets at the 
same time over a 25 arcmin diameter field of view, using fibres 
of 1.2" aperture. We observed a total of 15 hours on-source, sub- 
divided in 15 individual integrations of 1 hour duration. 



Fig. Q] shows a map of the observed region. We observed 
37 compact objects within 12' of NGC 1399 and with 18 < 
V < 21 mag (-13.4 < M v < -10.4 mag). This magnitude 
range covers the UCDs and overlaps the bright end of the glob- 
ular cluster luminosity function (Mieske et al. 2004). All targets 
have confirmed cluster membership from spectroscopic surveys 
(Drinkwater et al. 2000, Mieske et al. 2002 & 2004, Richtler et 
al. 2004 & 2008), except for the two objects closest to NGC 
1399, which were selected on the basis of their morphology 
from imaging from the ACS Fornax cluster survey (Jordan et 
al. 2007). 

For the observations we used the HR09A grism, which pro- 
vides an instrumental resolution of 8 km/s in terms of Gaussian 
cr (or 19 km/s in terms of FWHM) over a wavelength range 
5100 < A < 5400 A. This resolution allows us to reliably mea- 
sure velocity dispersions > 10 km/s. 

3. Data reduction 

3. 1 . Basic reduction 

In order to remove the instrumental signatures from the data, we 
used the publicly available GIRAFFE data reduction pipeline 
from ESQj. This pipeline performs bias subtraction, flat-field 
division, wavelength calibration, and spectrum extraction. As 
such, it creates a wavelength calibrated ID spectrum from a raw 
2D spectrum. The pixel scale in the wavelength calibrated ID 
spectrum is 0.05 A per pixel, slightly over-sampling the instru- 
mental scale of 0.08 A per pixel. The RMS of the wavelength 
solution was of the order 0. 1 A. The instrumental resolution — 
resulting from the four pixel FWHM of the fiber's spatial profile 
— in terms of Gaussian cr is ~0.14 A. This instrumental reso- 
lution corresponds to a velocity dispersion of <x~8 km/s in the 
5 100-5400 A wavelength regime. 

Given the multiplexing capability of FLAMES/GIRAFFE, 
we also obtained 23 sky spectra in each exposure by assigning 
unoccupied fibres to empty sky positions. These spectra were 
reduced identically to the science spectra. We combined the 23 
sky spectra in each exposure to 1 master sky spectrum, using the 
IRAF task scombine in the ONEDSPEC package. This com- 
bined sky spectrum was subtracted from each single, calibrated 
ID object spectrum. 

From this we obtained 15 sky-subtracted, calibrated ID 
spectra for each of the 37 compact objects observed. Those 
single spectra were corrected to heliocentric velocity using the 
IRAF tasks rvcorrect in the RV package and dopcor in the 
ONEDSPEC package. The velocity shift between the 15 indi- 
vidual spectra due to shifts in wavelength calibration was very 
small (< 2 km/s). We combined the 15 registered single spec- 
tra using the IRAF task scombine. For this we normalised the 
intensity of the spectra to their mode and applied a 3.5 <x aver- 
age sigma clipping algorithm. The resulting S/N per pixel in the 
combined object spectra ranged between 5 and 35. 

In order to have template spectra for measuring the inter- 
nal velocity dispersion, we also observed several dozen red gi- 
ant stars in the Milky Way globular cluster u> Centauri in a sin- 
gle FLAMES/GIRAFFE pointing. These stars are of late spec- 
tral type (typical temperature 5000 K), cover a metallicity range 
-2.3 < [Fe/H] < -0.5 dex and have magnitudes around V — 12 
mag (van Loon et al. 2007). We used the same instrument setting 
and reduction procedures as for the science targets. The internal 

1 http://www.eso.org/projects/dfs/dfs-shared/web/vlt/vlt-instrument- 
pipelines.html 
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Fig. 1. Left: Map of the observed region in the Fornax cluster. Hexagons indicate compact objects observed with FLAMES. Green 
squares indicate objects for which reliable velocity dispersions were measured. Large circles show targets with HST imaging that 
were successfully observed with FLAMES. Small dots are all known compact cluster members with V < 22 mag (My < -9.4 mag), 
large dots are those with My < -1 1 mag, the approximate magnitude division between UCDs and GCs. The dotted circles indicate 
Fornax cluster member galaxies from the Fornax Cluster Catalog (FCC; Ferguson 1989), for which the circle size gives the radius 
at which /iy = 25 mag / arcsec 2 . The asterisk marks the location of NGC 1399. Right: Magnitude distribution of observed sources. 
The solid histogram refers to those objects for which reliable velocity dispersions were measured. The dotted histogram refers to all 
sources included in the fibre configuration (see the left panel). 



line width of the giant stars is negligible compared to the in- 
strumental resolution. With 5 minutes on-source integration, we 
reached S/N ratios between 50 and 100. We used the 14 highest 
quality spectra as templates for the dispersion measurements. 

In Fig. [2] we show examples of two object spectra and one 
template spectrum. 

3.2. Dispersion measurements 

The internal velocity dispersion, <x, of each compact object was 
measured by cross-correlating its spectrum with various tem- 
plate spectra (IRAF task fxcor; Tonry & Davis 1979). For this 
measurement we excluded the wavelength region around the 
very strong Mgb lines (/i res tframe < 5210 A), since the measured 
width in this region proved to be systematically larger than in 
the rest of the spectra (see Fig. [3}. Such an increased width in 
the very deep a element absorption features is likely caused by 
saturation effects, and has been found previously by other au- 
thors (e.g., Rejkuba et al. 2007 and Hilker et al. 2007), who 
also excluded this region from their measurements. We used the 
wavelength region 5210 < /l res tframe < 5390 A, which includes 
the many prominent Fe absorption features around 5325 A (see 
Fig.©. 

Prior to cross-correlation, we continuum subtracted the spec- 
tra. For this, we adjusted the continuum fitting order individu- 
ally for each source such as to yield the lowest order that gives 
satisfactory results. The peak position of the cross-correlation 
gives the relative radial velocity between object and template. 
The width, cr pea k, of the cross-correlation peak (Fig. |4]i is the 
quadratic sum of the intrinsic object line width caused by ran- 



dom stellar motion plus twice the instrumental line width (equal 
to the template line width): i.e., cr 2 eak = cr 2 ^ + 2 x cr 2 ns . By 
cross-correlating the un-broadened and continuum subtracted 
templates against each other, we measured the template's intrin- 
sic line width cr; ns to be 9.7 km/s with a very small scatter of or- 
der 0.4 km/s. The intrinsic line width cr bj of the object spectrum 

is then calculated as: cr ob j = ^cr 2 eak - 2 x <x 2 ns . Note that the 

factor 2 in front of cr 2 ns is necessary because both the object and 
template spectrum are broadened by the instrumental resolution 
(Dubath et al. 1992). In Fig. @] we show the cross-correlation 
results from the two objects whose spectra are shown in Fig. [2] 

We performed tests with artificially broadened template 
spectra and different low frequency Fourier filter cutoffs to assess 
the accuracy of the fxcor task in measuring cr peak . We found 
a slightly non-linear relation between input spectral width and 
width measured by fxcor. The best agreement was found for a 
low-frequency cutoff of k=3 (see Fig. |5). We adopted this cut- 
off for the Fourier filtering, and applied a residual correction as 
a linear function of <x bj — as indicated in Fig. [5] — to the mea- 
sured width of the science spectra. The residual correction is in- 
dependent of the S/N in the object spectra, which we tested by 
artificially degrading the broadened template spectra to a range 
of S/N values between 5 and 35 per pixel, representative for our 
compact object sample. From the tests with the template spectra 
we also found that the background value in the cross-correlation 
peak fit needs to be kept fixed at (see also Fig. Allowing 
the program to fit the background value led to consistently over- 
estimated widths. 
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Fig. 2. Continuum normalised spectra of one template and two 
science objects, shifted to the restframe. The wavelength region 
/irestframe < 521 OA containing the Mg features was excluded from 
the cross-correlation fit (see Fig.[3]and text). 



We accepted a reliable measurement of cr ob j for a given ob- 
ject if two conditions were met: (1) the average confidence level 
of the cross-correlation peak was R > 4; and (2) none of the 
template cross-correlations yielded an outlier in the template- 
object relative velocity. The first condition removed 9 sources 
from the main sample of 37 objects while the second condi- 
tion removed five more sources. Fig. Q] shows that the rejected 
sources are mostly close to the faint magnitude limit of our sur- 
vey. We note that the two brightest sources with unreliable mea- 
surements (see also Fig. [TJ are those that had been selected as 
UCD candidates based only on morphology from ACS imag- 
ing (ACS Fornax cluster survey, see Jordan et al. 2007). Both 
sources are located within 2' to the center of NGC 1399 and are 
the only objects in our target sample whose coordinates could 
not be tied to the USNO B2.0 system. We attribute their low flux 
level to an offset in relative coordinates with respect to the rest 
of our sample. 

A final sample of 23 reliable measurements is obtained, 
of which 15 have archival HST imaging available. This is the 
largest homogeneous set of UCDs for which dynamical masses 
have been derived. The resulting range of intrinsic velocity dis- 
persions is 9 < <r ob j < 36 km/s, with a mean of 24 km/s. These 
values are listed in Table Q] for the 15 sources with HST imag- 
ing, and in Table [3] for the 8 sources without HST imaging. A 
map and the magnitude distribution of the investigated compact 
objects is shown in Fig.Q] 

3.3. Mass modelling taking into account aperture effects 

To estimate the masses of the UCDs, we used the mass mod- 
elling algorithm outlined in Hilker et al. (2007). This includes a 
correction of the measured dispersion cr ob j to the true global dis- 
persion cr bj >cor , due to the fact that our measurements miss the 
contributions from stars outside the fibre aperture. 
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Fig. 3. The x-axis shows the Gaussian width cr Mg in km/s of 
the cross-correlation peak derived in the Mg region (5170 < 
^restframe < 5210 A). The y-axis shows the difference between 
this width and cr derived in the rest of the spectrum (5210 < 
^restframe < 5390 A). For cr Mg > 30 km/s, the cross correlation 
width in the Mg region is significantly broader than in the rest of 
the spectrum. 
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Fig. 5. Plot illustrating the tests performed with cross-correlating 
artificially broadened spectra. The input width, cruy, is plotted 
against the width, <t ut, measured using the fxcor task. The 
solid line shows the identity relation. The dotted line is a linear 
fit to (Tout as a function of <x IN . The measured values of cr are 
corrected for according to this relation. The values shown here 
are for the case of k=3 as low frequency cut-off. For k > 3, the 
deviations from a linear relation were larger. 
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Fig. 4. Plot showing the Fourier cross-correlation peaks for objects F-19 (UCD3, left) and F-22 (UCD48, right), whose spectra are 
shown in Fig. [2] The lower panels are magnified views of the overall cross-correlation results from the upper panels. In the lower 
panels, the Gaussian fit to the cross-correlation peak is indicated by a dotted line. See text for more details. 



Table 2. Columns 2 to 6 of this table show the structural parameters of the King profile fits for the 15 UCDs from Table Q] (see 

For the brightest 



also Evstigneeva et al. 2008). The King profiles are parameterized as follows: Hf) - Ir, ■ > ■ r 

source (F-19), a composite King+Sersic profile was necessary to provide a satisfactory fit. The Sersic profile is parametrized as 



follows: I(r) = I e g exp 



The parameters of the Sersic profile are indicated in columns 7 to 9. 



IDflames 


/io [mag/arcsec 2 ] 


t c [pc] 


h [pc] 


c 


a 


Heg [mag/arcsec 2 ] 


r eff [pc] n 


F-19 


16.03 


4.90 


230.6 


1.67 


2 


21.34 


118.9 1 


F-24 


15.11 


3.03 


4501.2 


3.17 


3.32 






F-l 


14.81 


2.23 


487.1 


2.34 


1.23 






F-5 




1.22 


77.1 


1.80 


2 






F-12 


13.66 


1.24 


72.9 


1.77 


1.25 






F-ll 




1.70 


0.98 


49.1 


2 






F-9 




1.48 


3.14 


95.0 


2 






F-17 




1.70 


0.90 


45.0 


2 






F-7 


16.23 


7.03 


96.3 


1.14 


2.79 






F-22 


11.67 


0.39 


102.8 


2.52 


1.20 






F-6 




1.57 


2.29 


85.0 


2 






F-34 




1.51 


1.30 


42.2 


2 






F-51 




2.45 


0.49 


138.8 


2 






F-53 




1.95 


0.91 


81.5 


2 






F-59 




1.48 


1.97 


59.5 


2 







The mass modelling involved the following steps: 

1. The observed, PSF-deconvolved luminosity profile from 
HST photometry was parameterized by the best-fitting den- 
sity law. For most UCDs a satisfactory fit was achieved 
with a King or generalized King profile. Only UCD3 
(IDflames =F- 19) required a two component King+Sersic 
function to be fitted well (see also Evstigneeva et al. 2007). 



The profile parameters for all 15 sources are shown in 
Table |2] 

2. The 2-dimensional surface density profile was deprojected 
by means of Abel's integral equation into a 3-dimensional 
density profile. 

3. The cumulated mass function M(< r), the potential energy 
<p{r) and the energy distribution function /(£) were calcu- 
lated from the 3-dimensional density profile. 
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Table 3. Measured velocity dispersions cr bj for the Fornax UCDs without available HST imaging. IDs are from literature sources 
as in Table[TJ *This source designation is from Richtler et al. (2008). 



IDflames 


ID, , , 

literature 


RA [2000] 


DEC [2000] 


V [mag] 


Mv [mag] 


v,-ad.O [km/s] 


cr ob j [km/s] 


F-3 


UCD27 DW 


3:38:10.4 


-35:24:06.2 


19.7 


-11.7 


1626 


31.3(1.5) 


F-18 


UCD44 DW 


3:38:42.0 


-35:33:13.0 


19.7 


-11.7 


2024 


19.1 (1.4) 


F-23 


UCD49 D w 


3:39:20.5 


-35:19:14.2 


19.7 


-11.7 


1480 


21.9(1.4) 


F-2 


FCOS 2-2153 


3:38:06.5 


-35:23:04.0 


20.0 


-11.4 


1426 


18.7 (1.4) 


F-8 


FCOS 0-2066 


3:38:23.2 


-35:20:00.7 


20.1 


-11.3 


1414 


25.7 (1.4) 


F-40 


92.099* 


3:37:52.5 


-35:28:57.9 


20.7 


-10.7 


1497 


27.3 (1.4) 


F-60 


FCOS 2-2100 


3:38:00.2 


-35:30:08.2 


20.9 


-10.5 


871 


24.3(1.9) 


F-64 


FCOS 1-2080 


3:38:41.4 


-35:28:46.6 


21.0 


-10.4 


1728 


24.7 (1.7) 



4. Finally, an N-body representation of the UCD was created by 
using the deprojected density profile and the energy distribu- 
tion function. For every model, 100.000 test particles were 
distributed and their x, y and z positions and corresponding 
v x , v y and v z velocities were given as output. 

After generating the UCD model, the velocity dispersion as 
seen by an observer was simulated. In doing so, the following 
steps were performed: 

1 . All test particles are convolved with a Gaussian whose full- 
width at half-maximum (FWHM) corresponds to the ob- 
served seeing. 

2. The fraction of the 'light' (Gaussian) that falls into the fibre 
aperture (1'.'2 for FLAMES) at the projected distance of the 
observed object (19 Mpc) is calculated. 

3. These fractions are used as weighting factors for the veloc- 
ities. All weighted velocities that fall into the fibre region 
are then used to calculate the 'mimicked' observed velocity 
dispersion <r mo d- 

Iteratively, the total 'true' mass of the modelled object, M true , 
that corresponds to the observed velocity dispersion, <x t, s is then 
determined by scaling a first 'guess' mass, M guess , via the for- 
mula M true = M guess ' (o"obs/o" mo d) 2 - 

The masses, mass-to-light ratios, global and central velocity 
dispersions and the main model parameters derived in this way 
are listed in Table [TJ Note that the three brightest UCDs in our 
study were also observed by Hilker et al. (2007) with UVES. 
The M/L ratios derived for them in that study are indicated in 
Table [TJ The error ranges of the three estimates overlap, and the 
average ratio between our M/L values and those from Hilker et 
al. (2007) is 0.90 ± 0.16, consistent with unity. 

4. Results 

With the total mass derived from the dynamical modelling, we 
calculated the optical mass-to-light ratio MfLy, using the V-band 
photometry from the wide-field imaging data presented in Hilker 
et al. (2003) and Mieske et al. (2006, 2007a). Based upon the 
modelling algorithm of Hilker et al. (2007), we also calculated 
the central velocity dispersion <x for all sources with HST imag- 
ing available. The ratio ^ was 1 .23 on average with a scatter of 
0.07. For those eight sources with reliable cr measurements but 
without HST data (Table [3]), we assumed an average correction 
factor of 1.23 to include them in an analysis of their locus in the 
My - ctq plane. This plot is shown in Fig. [6] 

In the My - cr plane, the relation for globular clusters and 
the extrapolation of the Faber-Jackson relation for luminous 
elliptical galaxies (Faber & Jackson 1976) intersect at about 



My =: -10 mag. In the luminosity regime of our sample of 
Fornax UCDs, both extrapolations bifurcate. We can therefore 
roughly subdivide our sample into objects closer to the extrap- 
olation of GC relation, and objects closer to the extrapolation 
of the Faber-Jackson relation (which also happens to match the 
compact elliptical galaxy M32; Evstigneeva et al. 2007). The 
three brightest UCDs (My < -12 mag) are clearly more con- 
sistent with the Faber-Jackson relation while fainter UCDs are 
preferentially closer to the GC relation. We find that the pro- 
jected clustercentric distance of objects that are more consistent 
with the Faber-Jackson relation is ~ 60 + 30% larger than that of 
sources more consistent with the GC relation. These findings are 
consistent with the hypothesis that our sample consists of both 
objects associated to the globular cluster system of NGC 1399, 
and objects with more complex dynamical formation history, be- 
ing associated more to the overall cluster potential. 

In Fig. [7] we plot metallicity Z/H against M/L ratio. The 
metallicity Z/H is derived directly for some sources from pre- 
vious spectroscopy (Mieske et al. 2006), for others derived from 
their (V-I) colour (Mieske et al. 2007a), using the calibration 
of Kissler-Patig et al. (1998). This calibration was shown to be 
accurate to within 0.1-0.2 dex (Mieske et al. 2006) for old stel- 
lar populations with [Fe/H] > -1.0 dex. In the plot we indicate 
SSP predictions from Bruzual & Chariot (2003) and Maraston et 
al. (2005) for M/L ratios of populations with solar [a'/Fe] abun- 
dances, with ages between 5 and 13 Gyrs. The former models 
assume a Chabrier IMF (Chabrier 2003), while the latter models 
assume a Kroupa IMF (Kroupa 2001). We note that both IMFs 
are very similar and do not account for the difference in predicted 
M/L ratios at fixed age and metallicity (see also Dabringhausen 
et al. 2008). Rather, it is the choice of different stellar evolution- 
ary codes which leads to the 20% differences between Bruzual 
& Chariot and Maraston M/L predictions. Most of the M/L data 
points are consistent to within their errors with the theoretical 
predictions assuming a canonical IMF, with three sources show- 
ing somewhat elevated M/Ls. 

In the two top panels of Fig. [8] we plot My and mass vs. the 
M/L ratio for the 15 sources with HST imaging. We indicate the 
faint magnitude limit of our survey, which translates into a mass 
dependent M/L sensitivity limit. To test whether the rise of M/L 
with mass generally observed in the regime 10 5 < M/Mq < 10 s 
can be traced by our data, we fit a linear relation to the distribu- 
tion of mass vs. M/L ratio. We find a slope different from at 
the 2.8cr level. The significance of the slope was calculated by 
random resampling of the data points around the fitted relation. 
For this re-sampling, the scatter of the data points around the 
relation was used, given that it was about 25% larger than the 
average error of the data points. 

There are two caveats regarding the interpretation of this 
~3cr slope. The first caveat is the mass dependent M/L sensitiv- 



Table 1. Table showing the measured velocity dispersions, cr bj, for the Fornax UCDs with available HST imaging, plus the true global and central velocity dispersions <x bj iC or 
and <t , derived from modeling the light distribution of these sources within the FLAMES fibre aperture (see text and Hilker et al. 2007 for further details). We also indicate the 
derived total mass, the optical mass-to-light ratio M/Ly, and the effective radius r e ff. Coordinates are taken from USNO. The first error of M/Ly corresponds to the measurement 
uncertainty including that of the apparent magnitude. The second one is the systematic error arising from an assumed 0.15 mag uncertainty in the distance modulus of (m- 
M)=31.4 mag to Fornax. The heliocentric radial velocities have precisions of ~3 km/s, and a global uncertainty of 1-2 km/s, tied to the heliocentric velocity of a> Centauri (233 
km/s, Dinescu et al. 1999) . The fourth-to-last column shows the dynamical M/L derived by Hilker et al. (2007) for the three brightest sources. The three last columns indicate 
the relaxation time t le iax. the acceleration parameter a, and the projected distance d to NGC 1399 (l'=5.5 kpc at the Fornax cluster distance). 'Designation from Drinkwater et al. 
(2003). "Literature ID refers to Mieske et al. (2004) for FCOS, and Firth et al. (2007) for the UCDxx DW sources. 
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Fig. 8. Top Left: Absolute magnitude of UCDs plotted against their M/L ratio in the V-band. The dotted line indicates the faint 
magnitude limit of our survey. Top Right: Masses of the UCDs from the left plot, plotted against their M/L ratio. The dotted line 
indicates the mass dependent upper limit of our M/L sensitivity, caused by our faint magnitude limit of My = -10.4 mag. There is a 
3<x correlation between the shown data points. Bottom left: Plot analogous to the top right panel. Now, the M/L ratio measurements 
have been normalised to the same (solar) metallicity (see text). For this, we assume the mean of the model predictions from Bruzual 
& Chariot (2003) and Maraston (2005) in Fig. [7] The horizontal lines indicate the M/L ratios from the model predictions for 13, 9, 
and 5 Gyrs (top to bottom). There is only a 1.5<x correlation between the data points. Sources marked by large squares are those 
closer to the Faber Jackson relation in Fig. [6] Bottom right: As in the plot on the left, but now the x-axis is relaxation time. 



5. Discussion 

5.1. M/L ratio measurements over a range of environments 

In order to further quantify the dependence of M/L ratio on mass, 
and to investigate M/L variation with environment, we combine 
our data for Fornax UCDs with M/L ratio measurements of other 
compact stellar systems from literature studies, see Fig. [12] This 



covers the regime of low-mass Milky Way globular clusters (M> 
5 x 10 4 M Q ) up to the most massive UCDs (M~10 8 M Q ). 

5.1.1. Revision of literature mass estimates 

In Fig.Q~2]we use revised dynamical mass estimates of the Cen A 
globular clusters from Rejkuba et al. (2007) and the DGTOs 
from Hasegan et al. (2005). Those estimates are obtained by ap- 
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Fig. 6. Absolute magnitude of UCDs plotted against their central 
velocity dispersion cr . The solid line is the extrapolation of the 
relation defined by Galactic GCs (McLaughlin & van der Marel 
2005), the dashed line is the extrapolation of the Faber-Jackson 
relation for luminous elliptical galaxies (Faber & Jackson 1976), 
which also fits the compact elliptical M32. Open symbols indi- 
cate objects without HST imaging, for which the average cor- 
rection factor between global cr and cr was assumed as derived 
from modelling of the 15 objects with HST imaging. 




Fig. 7. Z/H of the UCDs plotted against their M/L ratios. Solid 
(green) lines indicate stellar population models from Bruzual 
& Chariot (2003) for ages (from bottom to top) of 5, 9, and 
13 Gyrs, assuming a Chabrier IMF. Blue (dashed) curves are 
from Maraston et al. (2005) for the same age ranges, assuming 
a Kroupa IMF. The Z/H values are from Mieske et al. (2006), 
assuming a solar [a/Fe] abundance, as suggested in Mieske et 
al. (2007b). Red data points indicate sources for which spectro- 
scopic [Fe/H] estimates are available. For the remaining sources, 
[Fe/H] is estimated from their ( V-I) colour, using the transforma- 
tion of Kissler-Patig et al. (1998) (see also Mieske et al. 2006). 



plying the same mass modelling and aperture simulations as pre- 
sented in Sect. |3.3| and in Hilker et al. (2007). To this end, we as- 
sumed King models as representations of the light profiles with 
the projected half-light radii and King concentrations as given in 
the two papers. We include new measurements from five addi- 
tional Cen A globular clusters (Rejkuba et al. 2007), using struc- 
tural parameters provided by M. Gomez from observations with 
IMACS on Magellan (private communication). 

It turned out that the modelled masses (and thus mass-to- 
light ratios) of the Cen A globular clusters are up to 30% lower 
than estimated from the virial mass estimator, whereas most of 
the modelled DGTO masses are in general higher (up to 60%). 
Figure [10] compares our modelled masses with the dynamical 
estimates given in the literature. The reason for the discrepancies 
is the different treatment of the aperture corrections that have to 
be applied to the observed velocity dispersions. 

In case of the Cen A data the authors estimated the aper- 
ture corrections to be a few percent, but they preferred to as- 
sign these corrections only to the total error budget, and instead 
used directly the observed (cr obs ) values with the virial estima- 
tor to derive masses. However, at the distance of Cen A, 3.84 
Mpc (Rejkuba 2004), the ratio between the projected half-light 
diameter (4-48 pc) and the slit width (1" = 19 pc) is such that 
aperture corrections can not be neglected. A large fraction of 
the light of the most extended massive clusters lies outside the 
slit area which prohibits the measurement of a global velocity 
dispersion. The ratios cr tot /cr b s and cr /cr obs as function of pro- 
jected half-light radius are shown in Fig.fTTI 
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Fig. 9. Left panel: Gravitational acceleration, a = ^j^> plotted 

against normalised M/L ratio. None of the objects are in the low 
acceleration regime a < 1.2 * 10~ 10 m s~ 2 where MOND has 
been postulated to hold. Right panel: Projected distance to NGC 
1399 plotted against normalised M/L ratio. 



In case of the DGTO data, the authors derived the dynamical 
mass from the King mass estimator (e.g., Dubath & Grillmair 
1997), for which the core radius and the central velocity disper- 
sion have to be known. Hasegan et al. (2005) corrected the mea- 
sured velocity dispersion to the central one by "scaling upward 
to account for the blurring of the actual velocity dispersion pro- 
files within the ESI slit" for which the knowledge of the intrinsic 
light profile from HST imaging has to be known. The Virgo clus- 
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Fig. 12. Top left panel: Mass vs. M/L for the isolated compact objects of this paper (green dots) plus literature values of compact 
objects in CenA (Rejkuba et al. 2007, cyan), Virgo (Hasegan et al. 2005, red; Evstigneeva et al. 2007, magenta), Fornax (Hilker et al. 
2007, green asterisks), and Milky Way globular clusters (McLaughlin & van der Marel 2005, black dots). The vertical dashed line 
indicates the approximate mass where the relaxation time is equal to one Hubble time (see also the right panels). Top right panel: 
Like the left panel, but plotting relaxation time instead of mass on the x-axis. Bottom left panel: Like the top left panel, but here all 
M/L ratio estimates have been normalised to solar metallicity (see text). The horizontal dotted lines indicate the M/L ratios expected 
for single stellar populations of age 13, 9, and 5 Gyrs (from top to bottom) based on Bruzual & Chariot (2003) and Maraston et al. 
(2005). The small (black) dots indicate the present-day M/L ratios of the Galactic GCs if they would not have undergone dynamical 
evolution (see text). Bottom right panel: Like the left panel, but plotting relaxation time instead of mass on the x-axis. 



ter is far enough away, 16. 1 Mpc (Tonry et al. 2001), that most of 
the light of the DGTOs, which have half-light diameter between 
7 and 58 pc, falls into the slit (with width of 0775 - 58 pc). 
One tends to measure the global velocity dispersion for most of 
the DGTOs. Thus, the correction to central velocity dispersions 
should be larger for smaller objects, since for those the mea- 
sured dispersion is closest to the global one. This, however, is 



opposite to the trend of the corrections applied by Hasegan et al. 
(2005) (see their Tables 4 and 5), such that on average the masses 
derived by Hasegan et al. (2005) have to be corrected upwards 
(see Fig. [TOt. We will take our modelling results for the further 
discussions. The revised model masses will also be included in 
an upcoming paper (Hasegan et al. 2008, in preparation), along 
with M/L measurements of newly discovered Virgo DGTOs. 
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Fig. 10. The effect of proper aperture corrections on dynamical 
mass estimates. The ratio between modelled mass and dynami- 
cal mass as given in the literature is plotted vs. modelled mass. 
The masses of Cen A GCs from Rejkuba et al. (2007) (solid cir- 
cles) were previously overestimated by up to 30%, whereas the 
masses of DGTOs in Virgo from Hagegan et al. (2005) (open 
triangles) were mostly underestimated. The errorbars reflect the 
uncertainties given for the literature values. 
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Fig. 11. Corrections from observed to total and central velocity 
dispersions for Cen A GCs based on our mass modelling. The 
upper panel shows the ratio between total and observed velocity 
dispersion, while the lower panel shows the ratio between cen- 
tral and observed velocity dispersion. The dashed vertical line 
indicates half the width of the slit aperture. 



5.1 .2. A trend of M/L with mass and relaxation time 

In Fig. [12] we plot both, mass and relaxation time, vs. direct and 
normalised M/L ratios for our compiled sample of compact stel- 
lar systems. In Table [5] the M/L ratios, metallicities and sizes 
for these objects are shown. When fitting a linear relation to re- 
laxation time as a function of mass, we find that relaxation time 
equal to 1 Hubble time is reached at a mass of about 2x10 6 Mq. 
In the following we adopt this mass as an approximate limit be- 
tween globular clusters and UCDs (see also Mieske & Kroupa 
2008). 

Fig- [J~2]shows a clear rise of the M/L ratio for masses above 
^2 x 10 6 Mq. When correcting the M/L ratio measurements 
for their metallicity dependence, this is still clearly visible (see 
also Mieske & Kroupa 2008; Dabringhausen, Hilker & Kroupa 
2008), formally significant at the 8 <x limit. The M/L ratio distri- 
bution of objects below and above the 2 x 10 6 M0 limit stems 
from the same parent distribution with only the 2.7x 10 _1 1 proba- 
bility, according to a KS test. The average normalised M/L ratio 
for globular clusters is 2.70 + 0.17, while it is 5.44 + 0.37 for 
UCDs. The rise in M/L corresponds to a 40% increase of nor- 
malised M/L ratio per mass decade: dl ° A f ^ = 0.147 + 0.019 
(applying a 3cr clipping). A separation between UCDs and GCs 
at freiax =1 tHubbie leaves the mean M/L ratios of GCs and UCDs 
unchanged with respect to the mass cut at 2 x 10 6 Mq. For 
the separation at f re ] ax =1 t Hu bbie> the M/L distributions have a 
common parent distribution at a probability of 1.4xl0~ 10 . M/L 
scales with relaxation time almost in the same way as with mass: 



dlog(M/L)) 



= 0.150 ±0.021. 



d/ogOrclax) 

Apart from studying the relative difference in M/L ratio be- 
tween GCs and UCDs, it is also important to compare the M/L 
ratios with the model predictions on an absolute scale. While 
the average M/L of UCDs is ~40% above the 13 Gyr isochrone, 
M/L ratios of the Galactic globular clusters are below the 13 Gyr 
isochrone by the same factor. The mean age derived from their 
location with respect to the isochrones is 7-8 Gyrs, well below 
the typical globular cluster age of ~12-13 Gyrs, indicating that 
the input for the stellar population codes may not represent the 
globular cluster properties (see also Dabringhausen, Hilker & 
Kroupa 2008). In what follows, we therefore discuss the extent 
to which dynamical evolution may have changed M/L ratios for 
the compact stellar systems under investigation, and hence con- 
tributed to shaping the observed trend between mass and M/L. 



5.1 .3. Does dynamical evolution shape the M/L trend? 

Baumgardt & Makino (2003) showed that star clusters experi- 
ence a depletion in low-mass stars leading to a drop in M/L 
of up to 0.5 after about 0.8 dissolution timescales. This drop 
in M/L corresponds to about 30% of the mean measured M/L 
ratio of Galactic globular clusters. Dynamical evolution could 
hence be responsible for the lower M/L ratios of GCs, provided 
that their dissolution timescale is comparable to or smaller than 
a Hubble time. To estimate the dissolution timescale ?di ss for a 
typical Galactic globular cluster, we apply equation (6) from the 
recent study of Lamers et al. (2006), in which the photometric 
evolution of dissolving star clusters in the Galaxy's gravitational 
potential is investigated. 



t diss =6.60xlO z (^) 



2/ % 653 0.967-0.00825x/o£(M,/10 4 ) 



X t 



(4) 



M, is the initial cluster mass. The time-scale to depends on 
the tidal field of the environment. Lamers et al. (2006) adopt 
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Fig. 13. Normalised M/L ratio of the sample of compact stellar 
systems from Fig. Q~2]plotted against their effective half-light ra- 
dius ri, in pc. Small (red) dots indicate GCs, defined as compact 
stellar systems with M< 2x 1O 6 M0. Large (blue) dots are UCDs, 
defined as compact stellar systems with M> 2x 10 6 Mq. The ver- 
tical dashed line indicates the radius where the relaxation time is 
equal to one Hubble time. 



to =21.8 Myr, which is valid for a circular orbit in the Galaxy 
at 8.5 kpc radial distance (Baumgardt & Makino 2003). The 
median galactocentric distance of the Galactic GCs plotted in 
Fig. Q~2] is 9.2 kpc (McLaughlin & van der Marel 2005), very 
close to the assumed 8.5 kpc. Therefore, we also adopt to =21.8 
Myr for estimating fd; ss . Evaluating equation [4] we obtain 
fdiss -50 Gyrs for 10 5 M©, and f diss ^250 Gyrs for 10 6 M . That 
is, dynamical evolution should not have changed the primordial 
M/L ratios of compact stellar systems with masses above 10 6 
Mq, including the UCDs. For Galactic GCs in Fig.|4l masses are 
between 10 4 and 10 6 Mq, such that their dissolution timescales 
are closer to, and in some cases below, a Hubble time. Assuming 
that the GC M/L ratio decreases linearly up to a difference of 0.5 
after 0.8 dissolution timescales (Baumgardt & Makino 2003), 
we plot in the bottom left panel of Fig.Q~2]the expected M/L ra- 
tios of Galactic GCs if there was no dynamical evolution. The 
corrections are small — on average about 5%, which is neglibile 
in the context of this discussion. Provided that the absolute scale 
of dissolution times derived by Baumgardt & Makino (2003) is 
applicable to the Milky Way GCs included in this study, their 
M/L ratios should not have notably decreased due to dynamical 
effects from their initial value. 

The stellar population models used here indeed appear to over- 
estimate the M/L ratios of globular clusters with a canonical 
IMF by ~40%, and on average under-estimate the M/L ratios 
of UCDs by about the same amount. 

5.2. Environmental dependence of M/L ratios 

From Fig. Q~2]it is evident that M/L ratio measurements of the 
Fornax UCDs fit well into the general trend of M/L increasing 



with mass. However, it is also interesting to note that we have not 
found a Fornax UCD with such extraordinarily high M/L ratios 
as the three DGTOs from Hasegan et al. (2005). Fornax UCDs 
cover the same mass range as Virgo UCDs, but their average 
M/L ratio is only 0.61 ±0.11 that of the Virgo UCDs, or 0.71 + 
0.08 when excluding the Virgo UCD with the highest M/L ratio 
(S999). 

Such a M/L ratio difference may arise from age differences. 
If Fornax UCDs have luminosity weighted ages around 7 Gyrs 
and Virgo UCDs ages around a Hubble time, both populations 
would be equally inconsistent with M/L ratio predictions from 
stellar populations for their age, indicating the presence of (bary- 
onic or non-baryonic) dark matter (Mieske & Kroupa 2008, 
Dabringhausen, Hilker & Kroupa 2008). Derivation of spectro- 
scopic ages from line abundances for Virgo UCDs have shown 
that they are most consistent with old ages around a Hubble 
time (Evstigneeva et al. 2007). For Fornax, the situation is less 
clear. Mieske et al. (2006) find indications for intermediate ages 
in Fornax UCDs from relating Hp to metallicity sensitive line 
indices, but these data were not calibrated to the Lick system. 
Accurate age determinations for a comprehensive sample of 
Fornax UCDs by spectroscopy or multi-band photometry will 
allow to draw firmer conclusions in this context. 

The lower M/L ratios of Fornax UCDs — provided they are 
not explained by age differences — may support the assumption 
of different competing formation channels for UCDs (Mieske et 
al. 2006). In that context, the higher M/L ratios of Virgo UCDs 
could be interpreted as being due to dark matter. Goerdt et al. 
(2008) show that under certain conditions, remnants of tidally 
stripped dwarf galaxies can maintain a significant amount of 
dark matter (although see Bekki et al. 2003 for a different view). 
The lower M/L ratios in Fornax UCDs could then be explained 
by them being stellar super clusters formed without dark matter 
(Fellhauer & Kroupa 2002, Mieske et al. 2006). 

5.3. Fundamental scaling relations of stellar systems: how 
do UCDs and GCs fit in? 

Here we analyse how compact stellar systems such as UCDs 
and GCs fit into fundamental scaling relations for more extended 
stellar systems. In motivating this analysis, we show in Fig. [T3] 
the half-light radius of the compact stellar systems from Fig.fT2l 
plotted against their M/L ratio. We mark the radius at which the 
relaxation time is equal to one Hubble time, which is roughly 7 
pc (log(r/,)=0.82; or log(r/,/kpc)=-2.18). As also seen in Fig.fTSl 
this limit nicely marks the rise of M/L ratios between the regime 
of GCs and UCDs. 

In Fig.[14]we show the location of all compact stellar systems 
from Fig.[T2]in the so-called 'fundamental manifold' (Zaritsky 
et al. 2006a & 2006b, Zaritsky et al. 2008). The fundamental 
manifold concept aims at a unifying empirical description of the 
structural and kinematic properties of stellar spheroids. It relates 
the effective radius r e to velocity dispersion <x and effective I- 
band surface brightness I e . In their studies, Zaritsky et al. show 
that stellar spheroids from the scale of galaxy clusters (~ 10 5 pc) 
down to the scale of dwarf elliptical galaxies (~ 10 2 pc) appear 
to form a common sequence in this manifold. With the data from 
Fig.Q~2]we can extend these considerations down to the smallest 
stellar systems (~ 10° pc). 

To derive I e for the compact stellar systems, we use (V-I) 
measurements where available, and otherwise convert [Fe/H] to 
(V-I) using the calibration relation of Kissler-Patig et al. (1998). 
In the plots we adopt a simple mass limit of 2xl0 6 M to sepa- 
rate UCDs from GCs (see Fig. [HI. 
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Zaritsky et al. (2006a, 2006b, 2008) define the fundamental 
manifold relation as 

log(r c ) = 2 x log(cr) - log(I c ) - log(T^) (5) 

In this formulation, log(T^) is the effective mass-to-light 
ratio parametrized in terms of log(cr) and log(I e ). That is, 
log(Tg) = log(Tg(log(cr), log(I e )). The parametrization is de- 
termined from a fit of dynamically derived M/L ratio (using 
the virial theorem) as a function of log(cr) and log(I e ) (see 
Zaritsky et al. 2008). It is clear that the exact functional shape 
of log(Yg(log(cr), log(I e )) — and hence the location of the fun- 
damental manifold — depends on which stellar systems are in- 
cluded in the fit. For example, the original formulation of the 
manifold (Zaritsky et al. 2006a) does not include the heavily 
dark matter dominated Local Group dwarf spheroidal galaxies 
in the fit (Zaritsky et al. 2006b, Simon & Geha 2007). A revised 
formulation extending to the very large M/L of the dSphs was 
presented in Zaritsky et al. (2008; see Table 1 of that paper). 

We are interested in the link between UCDs and canonical 
galaxies on the one hand, and on the relation between UCDs and 
star clusters on the other hand. Therefore, we show in Fig. [14] 
the location of UCDs and GCs with respect to three different 
formulations of the fundamental manifold. The functional form 
of log( Yg(log(<r), log(I e )) for these three representations is indi- 
cated in Table [4] 

The first formulation does not include Local Group dwarf 
galaxies to the fit of log(Yg(log(<x), log(I e )), nor UCDs and GCs. 
This is the original manifold version from Zaritsky et al. (2006a). 
It is intriguing that in this formulation, UCDs extend the fun- 
damental manifold relation by more than a decade in r e , down 
to r f ~ 5-7 pc (note that r e is shown in units of kpc in 
Fig. PHI , Together with all other spheroids they follow a well 
defined linear function slightly inclined with respect to the orig- 
inal fundamental manifold, with a slope 0.92 + 0.01. The funda- 
mental manifold relation in this formulation breaks down only 
for the faintest dwarf spheroidal galaxies (My > —9 mag) 
and for globular clusters. A possible interpretation of this is 
that for the faintest dwarf spheroidals, the dark matter halo is 
de-coupled from the baryons (see also Zaritsky et al. 2006b) 
such that the continous relation of baryon packing efficiency vs. 
galaxy scale breaks down. Another possibility is that the faintest 
dwarfs are out of dynamical equilibrium. This aspect is closely 
related to the discussion of the origin of dwarf satellite galax- 
ies (dark-matter dominated cosmological substructure vs. tidal 
dwarf galaxy, see Kroupa et al. 2005). 

The second formulation includes all objects in Fig. [14] for 
the fitting log(Yg(log(<x), log(I e )) (Zaritsky, private communi- 
cation). Again, UCDs follow the manifold line, and only for 
log(r e ) < -2.2 they start to 'bend down'. Interestingly, at this 
radius also the transition between objects with relaxation times 
smaller and larger than a Hubble time occurs (Fig.[T3ll. Globular 
clusters show a large scatter, and do clearly not align along the 
manifold. 

The third formulation includes all objects for fitting except 
GCs and UCDs with log(r e ) < -2.2. This formulation hence 
excludes dynamically relaxed stellar systems from the fit. UCDs 
with log(r e ) > -2.2 align very well with the manifold, while 
globular clusters and smaller UCDs do not. 

Summarizing, UCDs with log(r e ) > -2.2 (r e >7 pc) appear 
to form a single family with larger stellar systems in the funda- 
mental manifold. The location of most GCs is inconsistent with 
the fundamental manifold extrapolated from larger stellar sys- 
tems. Dynamically un-relaxed stellar systems appear to form a 



single manifold, while the relaxed systems — due to their ad- 
vanced dynamical evolution — scatter very broadly around it. 



6. Summary and conclusions 

In this paper we have analysed the internal dynamics of 23 
ultra-compact dwarf galaxies in the Fornax cluster. The analy- 
sis is based on high-resolution spectroscopy obtained with the 
FLAMES spectrograph at the VLT. Our targets cover an approx- 
imate mass range of 10 6 < M < 10 8 M0 and a luminosity range 
-10.4 < My < -13.5 mag, overlapping the bright end of the 
globular cluster luminosity function. We also compare the dy- 
namical properties of UCDs and GCs, and put them into the con- 
text of fundamental scaling relations defined for larger galaxies. 
We obtain the following results: 

1. In the My-cr plane, we find that UCDs with My < -12 mag 
are consistent with the extrapolation of the Faber-Jackson 
relation for luminous elliptical galaxies. For My > -12 
mag, most objects are located closer to the extrapolation to 
brighter luminosities of the globular cluster My-cr relation. 

2. We derive dynamical M/L ratios for those 15 of the 23 UCDs 
for which HST archival imaging is available, taking into 
account aperture effects in the spectroscopy (Hilker et al. 

2007) . Three out of the 15 UCDs have dynamical M/L ra- 
tios too high to be explained by canonical stellar populations, 
but we do not find Fornax UCDs with M/L ratios as extreme 
as found for some Virgo UCDs (Hasegan et al. 2005). At a 
given metallicity, Fornax UCDs have on average 30 to 40% 
lower M/L ratios than Virgo UCDs. 

3. We normalise the dynamical M/L ratios of the 15 Fornax 
UCDs to solar metallicity, using predictions from stellar pop- 
ulation models (Bruzual & Chariot 2003, Maraston 2005). 
We find no significant correlation between normalised M/L 
ratio and mass or relaxation time for our Fornax UCD sam- 
ple. We do not find a dependence of normalised M/L ratio on 
projected clustercentric distance. 

4. We add our new measurements for 15 Fornax UCDs to the 
available data on M/L ratios of compact stellar systems in the 
broader mass range 10 4 < M < 10 8 Mq. We include Galactic 
globular clusters and UCDs in Virgo, CenA and Fornax. We 
re-analyse dynamical mass estimates of UCDs in the Virgo 
cluster (Hasegan et al. 2005) and the CenA group (Rejkuba et 
al. 2007), using our modelling algorithm (Hilker et al. 2007) 
to correct for aperture effects in the spectroscopy. The cor- 
rections for the global velocity dispersion cr are of the order 
of 5-10%. We also provide previously unpublished M and 
M/L estimates for 5 CenA compact objects. 

5. We find a clear break in the distribution of normalised M/L 
ratios at a characteristic mass of ^2xlO 6 M0, which roughly 
corresponds to a relaxation time of one Hubble time (see also 
Dabringhausen, Hilker & Kroupa 2008, Mieske & Kroupa 

2008) . Objects more massive than this limit have normalised 
M/L ratios twice as large as objects less massive than this 
limit. In this context we suggest to separate UCDs from GCs 
by a mass limit of M^ 2 x 1O 6 M0 (see also Hasegan et al. 
2005). On average, the M/L ratios of UCDs are 40% above 
the expectations for a 13 Gyr stellar population with canon- 
ical IMF, while for GCs they are 40% below these expecta- 
tions. We find that the M/L ratios estimates of GCs are proba- 
bly only weakly biased (~5%) by their dynamical evolution, 
indicating that stellar population models indeed over-predict 
M/L ratios for compact stellar systems like GCs. 
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Table 4. Fitting coefficients for log(Yg(log(cr), log(I e )) in the three fundamental manifold formulations from Fig. [14] The functional 
form is log(Yjj) = cl + c2 * log(cr) + c3 * log(I e ) + c4 * (log(I e )) 2 + c5 * (log(cr)) 2 + c6 * log(o-) * log(I e ). 



Manifold formulation 


cl 


c2 


c3 


c4 


c5 


c6 


1 


2.75 


-1.70 


-0.295 




0.63 




2 


1.8974 


0.1896 


-0.9699 


0.1095 


0.1193 


0.02893 


3 


2.2397 


-0.3006 


-0.8726 


0.1159 


0.2827 


-0.0337 



6. UCDs extend the 'Fundamental Manifold' in its original for- 
mulation (Zaritsky et al. 2006a) by more than a decade in r e 
down to r e ~ 5 -7 pc. In this formulation, neither the faintest 
dwarf spheroidals {My > -9 mag) nor GCs lie on the man- 
ifold. When using also GCs, UCDs and dSphs to define the 
shape of the FM, UCDs with r e > 7pc and dwarf spheroidals 
align along the manifold, while GCs and smaller UCDs do 
not. This characteristic scale of r e ^7pc also marks the tran- 
sition between compact stellar systems with relaxation times 
below and above a Hubble time. 

We suggest a defintion of UCDs as those compact stellar sys- 
tems with M>2xl0 6 M Q and 7<r e /pc<100. As such, UCDs are 
the smallest dynamically un-relaxed stellar systems. From their 
position in the 'Fundamental Manifold' they can be considered 
the small-scale end of the galaxy sequence. 

A key question about UCDs is whether they are of 'cosmo- 
logical' origin, hence related to compact low-mass dark matter 
halos. Their elevated M/L ratios can be interpreted as mark- 
ing the on-set of dark matter domination in small stellar sys- 
tems. However, dark matter can hardly be detected directly, such 
that observational efforts need to be directed towards verify- 
ing/excluding alternative scenarios, such as a variation of the 
IMF in UCDs (Mieske & Kroupa 2008). In parallel, theoreti- 
cal studies regarding the dynamical evolution of compact stellar 
systems embedded in dark matter halos are needed for the mass- 
size regime of UCDs. 

In this paper it has been found that Fornax UCDs have 30- 
40% lower dynamical M/L ratios than Virgo UCDs. A possi- 
ble explanation for this is that only Virgo UCDs have signif- 
icant fractions of dark matter. This may be explained by the 
dominance of different UCD formation channels in Virgo and 
Fornax (Mieske et al. 2006). A simple way to test the possibility 
of different dark matter fractions is to determine the luminos- 
ity weighted ages of Fornax and Virgo UCDs. Younger ages in 
Fornax UCDs of ~7 Gyrs would naturally explain the M/L ratio 
differences and imply similar dark matter fractions as in Virgo. 
Together with efforts to constrain the IMF shape in UCDs, such 
an observational study is the next logical step in understanding 
the puzzling nature of UCDs. 
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-2 2 

log(R e ) [kpc] 

Fig. 14. The location of the compact stellar systems from 
Figs. [12] and [13] with respect to three different formulations of 
the 'fundamental manifold' (Zaritsky et al. 2006a, Zaritsky et 
al. 2006b). Table [4] gives the fitting coefficients for the three 
formulations. Large blue dots are UCDs, small red dots are 
GCs. A global mass cut at 2xl0 6 M© is applied to separate 
UCDs from GCs. Cyan dots are large spheroids from the sam- 
ple of Zaritsky et al. (2006a), with open cyan circles indicat- 
ing faint Local Group dSphs with My > -9 mag (Zaritsky 
et al. 2006b, Simon & Geha 2007). Plot 1: This is the orig- 
inal formulation of the manifold (Zaritsky et al. 2006a), for 
which the fit of /og(Y{(log(cr), log(I e )) only includes large stel- 
lar spheroids with log{r e ) > -0.4, excluding the Local Group 
dSphs. Plot 2: The fit of fog(r{(log(cr), log(I e )) includes all ob- 
jects in the plot (Zaritsky private communication). Plot 3: The fit 
of logCf{{log(cr), log(I e )) includes all objects in the plot except 
UCDs with log(r e ) < -2.2 and GCs (Zaritsky, private communi- 
cation). See also Fig.[T3land text. The limiting log(r e ) is marked 
by a dashed vertical tick. 
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Table 5. This table gives masses, MfL v ratios, metallicities, M/L v ratios normalised to solar metallicity, half-light radii and velocity dispersions 
a for the compact stellar systems from Fig.[T2]to Fig. [14] The table is ordered by descending mass. Sources are this paper", Hilker et al. (2007)*, 
Rejkuba et al. (2007)'', Hasegan et al. (2005)'', Evstigneeva et al. (2007) c , McLaughlin & van der Marel (2005)' , Meylan et al. (2001) s , Barmby et 
al. (2007)'', de Marchi et al. (1999)'. For galactic sources (ID MW....; McLaughlin & van der Marel 2005), the quoted a measurements are close 
to central values. For extragalactic sources, cr is the global velocity dispersion, corrected for spectroscopy aperture losses (see text and Hilker et 
al. 2007). 
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Table 5. continued. 
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